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Abstract 

We investigate the behavior of entanglement entropy across the holographic 
p-wave superconductor phase transition in an Einstein- Yang-Mills theory with 
a negative cosmological constant. The holographic entanglement entropy is cal- 
culated for a strip geometry at AdS boundary. It is found that the entanglement 
entropy undergoes a dramatic change as we tune the ratio of the gravitational 
constant to the Yang-Mills coupling, and that the entanglement entropy does 
behave as the thermal entropy of the background black holes. That is, the 
entanglement entropy will show the feature of the second order or first order 
phase transition when the ratio is changed. It indicates that the entanglement 
entropy is a good probe to investigate the properties of the holographic phase 
transition. 



1 Introduction 

As a strong-weak duality, the AdS/CFT correspondence [1, 2, 3] provides a powerful 
method for studying a strongly interacting system through its gravity dual which is weakly 
coupled. Especially, it has been used widely to model basic phenomena in condensed mat- 
ter physics, such as superconductivity (superfluidity) [4], Nernst effect [5], and non-fermi 
liquid [6]. For more related studies, see, for example, Refs. [7, 8, 9] and references therein. 
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The physical picture behind the holographic superconductor model is as follows. As 
the simplest concrete model, consider an Einstein-Maxwell-scalar field theory with a neg- 
ative cosmological constant. At high enough temperature, the Reissner-Norstrom-AdS 
(RN-AdS) black hole with a trivial scalar field is stable. And the dual CFT is in a de- 
confined phase and describes a conductor phase. When one lowers the temperature of 
the black hole, the RN-AdS black hole becomes unstable, a new black hole solution with 
nontrivial scalar field is favored, which can describe a superconducting phase. The U(l) 
symmetry is spontaneously broken due to the nontrivial scalar field. The condensation 
of the scalar "hair" of the black hole gives a finite vacuum expectation value of the dual 
operator in the field theory side, which plays the role of order parameter in the holographic 
phase transition. The s-wave superconductor is described by the appearance of the scalar 
"hair" [4, 10, 11, 12], while the p-wave superconductor is characterized by the condensation 
of the vector "hair" [13, 14]. 

On the other hand, the entanglement entropy is expected to be a key quantity in our 
understanding some characterization of several aspects in many-body physics (see, for ex- 
ample, Refs. [15, 16]). For a given system, the entanglement entropy of subsystem with its 
complement is defined as the von Neumann entropy. In the spirit of AdS / CFT correspon- 
dence, a geometric proposal to compute the entanglement entropy has been presented in 
Ref. [17]. More precisely, consider a subsystem A of the total boundary system, the entan- 
glement entropy of subsystem A with its complement is given by looking for the minimal 
area surface 7^4 extended into the bulk with the same boundary dA of A (see Refs. [18, 19] 
for reviews) 

Area(7^) 

where Gn is the Newton's constant in the bulk. While various aspects of different holo- 
graphic superconductor models have been intensively studied (see, for example, Refs. [20, 
21, 22, 23, 24, 25, 26, 27, 28, 30, 31]), the study of entanglement entropy in the holographic 
phase transition is just in the early stage. Ref. [32] studied the behavior of entanglement 
entropy in a holographic s-wave superconductor model, while Ref. [33] discussed the case 
in the holographic insulator/superconductor phase transition. 

Note that a holographic p-wave superconductor (superfluid) with fully back reaction 
in the Einstein- Yang-Mills theory was constructed in Ref. [34]. This model is interesting 
not only because it is a holographic model to describe a p-wave superconducting phase 
transition, but also it contains a rich phase structure. There is a parameter a, the ratio 
of the gravitational constant to the Yang-Mills coupling, in this model. The p-wave su- 
perconductor phase transition is second order for small a, while it will become first order 
as a increases beyond a critical value. Therefore it is quite interesting to see the behavior 
of entanglement entropy in this model, in particular, to see how the entanglement entropy 
changes when the order of the phase transition changes. 

The aim of this paper is just to investigate the behavior of entanglement entropy in 
the holographic p-wave superconductor at finite temperature. The entanglement entropy 
is calculated for a straight strip geometry at AdS boundary by using of the holographic 
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proposal (1). We find that the behavior of entanglement entropy changes dramatically 
when the order of the phase transition changes. When the strip width is very large, i.e., 
7.4 probes deeply, the entanglement entropy is extensive as the thermal entropy of the bulk 
black hole, while in the opposite limit, the behavior perfectly fits the general form obtained 
from four- dimensional conformal theories. For the case with an intermediate strip width, 
by comparing the entanglement entropy and the thermal entropy of the bulk black holes 
during the whole process of phase transition, we see that they show the same behavior. 
This is an interesting and nontrivial result. As a result it shows that the entanglement 
entropy is a good probe to the holographic phase transition and that its behavior can 
indicate the appearance as well as the order of phase transition. 

The paper is organized as follows. In Section (2), we briefly review the holographic 
p-wave superconductor model and give the complete equations of motion to be solved. In 
Section (3), the fully back-reacted system is solved by shooting method and basic behaviors 
in equilibrium are described. In Section (4), we explore the behaviors of the entanglement 
entropy in the p-wave superconductor phase transition. The conclusion and discussions 
are included in Section (5). 



2 Gravity Background 

We begin with the Einstein- Yang-Mills theory in five-dimensional asymptotically AdS 
spacetime 

/I 19 1 

where k is the five dimensional gravitational constant connected with Gn by the relation 
2k 2 = IQttGn, g is the Yang-Mills coupling constant and L is the AdS radius. The SU(2) 
Yang-Mills field strength is 

F a , v = d,Al - d v Al + e^Al, (3) 

where v — (t,r,x,y, z) denote the indices of spacetime and a,b,c = (1,2,3) are the 
indices of the SU(2) group generators r a = a a /2i (a a are Pauli matrices). e abc is the totally 
antisymmetric tensor with e 123 = +1. The gauge field is given by A = A^^dx^ 1 . Here we 
define a parameter a = n/g which measures the strength of the back reaction. 

Following Refs.[14, 34], our ansatz for the metric and Yang-Mills field are chosen by 

ds 2 = -N(r)a(r) 2 dt 2 + -^-dr 2 + r 2 f(r)- A dx 2 + r 2 f(r) 2 (dy 2 + dz 2 ), (4) 

N(r) 

A = (p(r)r 3 dt + w{r)r l dx. (5) 
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The independent equations of motion in terms of the above ansatz are deduced as follows 

J 3r 2 N 2 a 2 3r 2 1 \r f N a 

= 2Ar -0 ), 
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where m(r) = ^p- — ^N(r) and " ' " denotes the derivative with respect to r. The event 
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We should demand 4>{th) = to have a finite form for gauge field at horizon. The 
asymptotical behavior of these fields near the horizon are 
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All coefficients in above expansions are constants and are related by the equations of 
motion (6). After substituting the expansion into (6), we find only five independent pa- 
rameters, i.e., {r H ,4>H \ w h \°h \jh }■ 

The ultraviolet (UV) asymptotic expansion near the boundary r — > oo behaves as 

X AO) , (0) . % 

_,(2) _(4) f (4) ^ 

(0) , m B (0) . O b (o) 

m = m B + + • • • > a = °b + — + •••' / = A + — + •••• 

To recover the pure AdS boundary, we need the boundary conditions = 1 and = 1. 
4>£p is the chemical potential /i and iw^ 1 is the source of the operator Jf . To spontaneously 
break the U(l) gauge symmetry and rotational symmetry, we should impose w^P = 0. 

After imposing boundary conditions, the equations of motion can be solved numerically 
via tuning the four independent parameters {4>h\ w h\ a H' '> fn°^} t° search for solutions that 
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meet the requirements a$ = 1, = l,w^ = 0. Notice that the above equations of 
motion (6) have four useful scaling symmetries [34] 

a ->• Act, -> A0, (9) 

f^Xf, co^\- 2 u, (10) 
r — > Ar, m — > X 2 m, L XL, — > A _1 0, a — >• Ace, (11) 
r — )• Ar, {t, x, y, 2;} — >• A~ 1 {t, x, y, z}, m — > A 4 m, — >• Acj, — >• A0. (12) 

Taking advantage of the scaling symmetries (9) and (10), we will first choose = 1, = 
1 in our shooting method, then use the two scaling symmetries again to set cr^ = 1, = 
1. The last two scaling symmetries allow us to set L — Th — 1- 

3 Thermodynamics and Phase Transition 

From the discussion in Section 2, for given a}, we can solve the equations of 

motion (6) by choosing (jy^ as a shooting parameter. After solving the coupled equations, 
we can obtain the condensate (Jf ) = ^rWg\ chemical potential /1 and total charge density 

p = ^r0g' ) by just reading off the coefficients w^j , (j)^ and 0^ from (8) respectively. 
However, there is an analytic black hole solution of (6) for vanishing u(r), which is just 
the RN-AdS black hole 

2 

0(r)=/x(l-^f), u(r) = 0, <j( r ) = f{r) = l, 

v 2 y 2 r H 1 - 2 - 2 
r 2 ' 3 r 2 ' 



, T . . 2 . a 2 u 2 rjr 1 a 2 u 2 9 1 4 . 9 
N(r) = --( + -^-r 2 H + -r|) + r 2 . 



This so called RN-AdS solution has vanishing u(r), thus corresponds to the normal phase. 
From AdS/CFT correspondence, the boundary thermal equilibrium states are dual to 
black hole geometries in bulk, and the Hawking temperature of black hole is considered as 
the temperature of the boundary theory [35]. From the metric ansatz (4), the Hawking 
temperature of the black hole is 



aN' 

T 



47T 



r=r H 
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The Bekenstein-Hawking entropy of this black hole is 



AG N h? 

where Ah denotes the area of the horizon and V = J dxdydz. 



St = = ^Vr H , (15) 
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We will work in the grand canonical ensemble where the chemical potential /i at the 
boundary is fixed. It is convenient to express physical quantities in scale invariant way. 
Under the scaling symmetry (12), the relevant quantities scale as follows 



fi A/i, T AT, p -> A 3 p, (Ji) -> A 3 (Jf). 



(16) 



Therefore, we choose the following scale invariant combinations to examine physics in the 
grand canonical ensemble 

T p (Jf) 

,3' 



(17) 



In gauge/gravity duality the grand potential Q of the boundary thermal state is identi- 




w(l/r) 



m(l/r) 




Figure 1: The configurations of gauge fields (<f), u) and metric functions (/, m) as a function 
of the inverse holographic coordinate 1/r for a = 0.447 at T/fi ~ 0.0215. 

fied with T times the on-shell bulk action in Euclidean signature. The Euclidean action 
must include the Gibbons-Hawking boundary term for a well-defined Dirichlet variational 
principle and further a surface counterterm for removing divergence 



S 



Euclidean 



J ^^^T^-W^^ (18) 



where h is the induced metric on the boundary r — > oo, and K is the trace of the extrinsic 
curvature. 
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This model has been numerically solved in Ref. [34] and shown that the order of the 
phase transition relies on the value of a. The transition is second order as a is less than 
a c = 0.365 ± 0.001, while it is first order for larger values than a c . We will re-solve the 
equations of motion (6) for completeness and for further discussion. Typical solutions for 
the metric and gauge field configurations are presented in Figure. (1), which are needed to 
calculate entanglement entropy in the next section. 

To compare the differences between the second order transition and first order transi- 
tion, we choose a = 0.316 < a c and a = 0.447 > a c as concrete examples. The conden- 




Figure 2: The condensation of vector operator Jf as a function of temperature for a = 0.316 
(left plot) and a = 0.447 (right plot). The vertical dashed line on the right plot represents 
the transition temperature. The right plot has a region where (Jf ) is multi-valued. 

sation of vector operator Jf as a function of temperature for a = 0.316 and a = 0.447 is 
displayed in Figure. (2). For the case a = 0.316, the condensate (Jf) appears at a particu- 
lar temperature T c ~ 0.0458/i. As the temperature is lowered, (Jf ) increases continuously. 
The critical behavior near T c is found to be ( Jf ) oc (1 — T/T c ) 2, which is the typical result 
from the mean- field theory. While for the case a = 0.447, we can see from Figure. (2) that 
the curve of the condensate has two branches when 0.905 < T/T c < 1.06. Therefore, the 
value of condensate has a jump at critical temperature T c ~ 0.0218/i, which represents a 
first order phase transition. 

To distinguish which branch is physical, i.e., thermodynamically favored, we need to cal- 
culate the grand potential Q. The values of grand potential Q are exhibited in Figure. (3). 
The RN-AdS solutions always exist for all temperatures, but it is only thermodynamically 
favored at high temperatures T > T c . As the temperature is lowered below critical value 
for each a, the grand potential Q from superconductor solution is smaller than the one 
from the RN-AdS solution, thus the superconductor phase is physically compared to the 
normal phase (RN-AdS solution). Therefore, there is a phase transition occurring at T c . 
For a = 0.316 case, the transition is second order. However, for a = 0.447 case, there is a 
characteristic "swallowtail" shape of the grand potential, signaling a first order transition. 

According to the two-fluid model, the total charge density p can be divided into two 
components p = p n +p s , where p n is the normal component, while p s is the superconducting 
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Figure 3: The grand potential Q as a function of temperature for a = 0.316 (left plot) and 
a = 0.447 (right plot). Trace the physical curve by choosing the lowest grand potential 
at a fixed T. The critical temperature T c is the point at which the superconductor phase 
begins to be thermodynamically preferred. In both plots, dashed blue curves are for the 
RN-AdS solutions, while the solid curves are for the superconductor solutions. 



component. In the holographic setup, the normal charge density p n is proportional to the 
t 3 part of the electric field at the horizon [14], which is given by in our units. 

Therefore the superconducting charge density is p s = p — p n , where p = "^t<$b * s the total 
charge density. We draw ^ and ^ versus temperature in Figure. (4). ^ and ^ as two 
functions of temperature are reminiscent of the temperature dependence of the superfluid 
and normal components of liquid He II as measured from in the torsional oscillation disk 
stack experiment. p s seems to vanish near T c like a power law, i.e., (T c — T) u . However, 
we find — goes to zero linearly here, while the experiment gives a value of v ~ 0.67. 



4 Entanglement Entropy 

After solving the equations of motion, we are now ready to calculate the entanglement 
entropy in this holographic model. Because of the arbitrary choice of the subsystem A, we 
can define infinite entanglement entropies correspondingly. However, we are here interested 
in a belt geometry with a finite width £ along the x direction and infinitely extending in y 
and z directions. 

To deal with the UV divergence, we assume that the subsystem A sites on the slice 
r — - where e — > is the UV cutoff. More specifically, •ja starts from x = — | at r = -, 
extends into the bulk until it reaches the minimum r = r*, then returns back to the AdS 
boundary r = - at x = +|. According to the proposal (1), we need to minimize the 
following area functional 

Area M = V, fj iJ^^^^ (19) 



s 



0.6 0.7 0.8 0.9 1.0 1.1 
T 



Figure 4: ^ and ^ are plotted versus temperature in a = 0.316 case. In our numerical 
calculation, we find — goes to zero linearly near T c , which is the same as the critical 
behavior obtained in the probe limit. 

where V2 — J dydz. The integrand can be considered as the Lagrangian with x direc- 
tion thought of as time. As the Lagrangian does not explicitly depend on "time x" , the 
Hamiltonian is conserved. Thus we can easily deduce the "equation of motion" that gives 
minimal area from (19) 

dr ryiVTJ 
dx f(r) 

where we demand that the surface is smooth at the turning point r = r*, i.e., dr/dx\ r = 
0. Integrated once, the belt width I can be fixed as 




(20) 



dr-=Pdr^- S= (21) 

2 J r . dr J rt ry/N(r)y/ r6 - r * 

Substituting (20) into (19), we finally obtain the entanglement entropy 

S E = ^ [■ 1 = % VA \ + Se) , (22) 

where the UV cutoff 1/e has been taken into consideration. The first term indicates UV 
divergent (e — > 0) and represents the "area law" [17, 36]. It can be deduced by plugging the 
UV asymptotic expansion (8) into (22). While the second term is independent of the cutoff 
and is finite, so this term is physical important. Following the discussion in Section (3), I 
and Se under the transformation (12) scale as I — > A _1 £, Se —> A 2 ^, so we will introduce 
the scale invariants 

I*, %• (23) 
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/// 

Figure 5: The entanglement entropy in a = 0.316 CclSG ; clS Sb function of belt 
width at fixed temperature. The curves from top to bottom correspond to T ~ 
0.0458/i, 0.0448/i, 0.0424/x, 0.0387/i, and 0.0294/1, respectively. 

We first focus on the case with second order phase transition, i.e., a < a c . In Figure(5) 
we plot the behavior of the universal part of the entanglement entropy Se as a function of 
belt width i by fixing the temperature. The curve at the top is at the transition temperature 
T c , which is identical with the RN-AdS case. We observe that the slope of the curve 
decreases as the temperature is lowered in superconductor situation as is expected that the 
lower the temperature is, the more the degrees of freedom will condense. This phenomenon 
can be seen much more clearly in Figure(6), which shows how the entanglement entropy 
evolves with temperature by fixing the belt width. We see that although the entanglement 
entropy is continuous at critical temperature T c , there is a discontinuous change in its slope 
at T c . This discontinuity may signal a significant reorganization of the degrees of freedom 
of the system, since some kind of new degrees of freedom, like the Cooper pair, would 
emerge in the new phase. 

We show the results for the first order transition case in Figure(7). The behavior of 
entanglement entropy at fixed temperature is quite similar to the second order transi- 
tion case. Se changes monotonously with respect to the belt width. This behavior is 
quite different from the result in Ref. [32], where a swallowtail shape appears in a region 
with finite belt width. It can be seen from the left plot in Figure (7) that the belt width 
monotonously decreases as the turning point r* increases. Actually, £ diverges logarithmi- 
cally as r* — > th at nonvanishing temperature as observed from (21). In fact in Ref. [32], 
the non-monotonousness of I with respect to the turning point r* is indispensable for the 
emergence of the kink. On this we will have more discussions in the last section of the 
paper. 

The behavior of Se as a function of temperature for fixed belt width is presented 
in Figure(8). Comparing with Figure(6), we can find a dramatic change in the a > a c 
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Figure 6: The entanglement entropy in a = 0.316 function of temperature at 

fixed belt width (left plot for l\x = 4 and right plot for 1\l — 6). The dashed blue curves 
are from the RN-AdS solutions, while the solid curves are from superconductor solutions. 
The physical curve is determined by selecting the lower entropy at a given T. 



case. The blue curve from normal phase is physical as T > T c , while the curve with the 
lowest entropy at a given temperature is preferred below T c . Therefore, there is an obvious 
jump in Se as well as its slope at critical temperature. It seems reasonable to expect an 
abrupt reduction in the number of degrees of freedom at T c since the condensate has a 
sharp jump at the critical point. 

In both two kinds of phase transition, we observe from Figure. (5) and Figure. (7) that 
Se exhibits linear behavior with respect to £ for large £. Indeed, we can find that in large 
£ji ~ iT limit, the main contribution of the integrals (21) and (22) to Se comes from the 
region near r = r* ~ r#. In addition to N(r) oc (r — rjj) near horizon at nonvanishing 
temperature, we can deduce the linear relation ~ i/i. Thus we obtain the entanglement 
entropy in this limit 

— id s V 2 £ 

K K 

where we have subtracted the UV divergent term. This equation is similar to the Bekenstein- 
Hawking entropy (15) and may seem to be surprised at first glance since the entropy (24) 
is proportional to the area of the subsystem A as opposed to the area law (1). A geometric 
interpretation made in Ref. [17] is that 7^4 will wrap a part of the black hole horizon as 
the increase of belt width and therefore is equal to the fraction of black hole entropy . 

In contrast, from Figure. (5), and Figure. (7), the value of seems to be power-law 
divergent as \il vanishes. In fact, the relationship between entanglement entropy and belt 
width in small l\i ~ £T limit is found to perfectly agree with a universal function in our 
numerical calculation 

& = ^(I -2^5?) = 1^-2^5). (25) 



S E ~ — V 3 V 2 £ ~ -TV, (24) 
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2 4 6 8 10 12 14 

Figure 7: The left plot shows the behavior of belt width as a function of turning point 
r* for a = 0.447 at T//J, ~ 0.0215. While the entanglement entropy as a function of belt 
width is presented in the right plot. 

Notice that entanglement entropy in four-dimensional conformal theories for a belt config- 
uration takes the universal form [37, 38] 

S A = Cj~ (26) 

where ( and C are numerical constants which depend on the details of a theory under 
consideration. Our IT — > result (25) has the same form as (26), which is reasonable since 
our metric solution is asymptotically AdS and 7.4 with small belt length can only probe 
the bulk sufficiently near the boundary r — > 00. As discussed in [39], there exist some 
universal crossover functions connecting the universal parts of the entanglement entropy 
to the thermal entropy. We try to construct the crossover function here as Se{£,T) = 
T n F(£T). In our situation, as £T — > 00, we obtain that n=2 and F(£ — > 00) ~ £, thus 
Se{£,T) behaves as the extensive thermal entropy. While as IT — > 0, we find n = 2 and 
F(£ — > 0) ~ ^. For the intermediate scale, the precise form of F(£) can only be obtained 
numerically. 

To get further understanding of the connection between the entanglement entropy and 
thermal entropy, it would be irradiative and instructive to compare the behaviors of them 
during the process of phase transition. Figure. (9) shows the thermal entropy St as a func- 
tion of temperature. The quite qualitatively similarity in Figure. (6), Figure. (8) and Fig- 
ure.^) is impressive and striking. Note that £fi = 4 and £/i = 6 in Figure(6) and Figure(8), 
this result is nontrivial since the belt width here is neither too large nor too small. The 
minimal surface 7,4 will not be very close to the horizon. Due to the lake of numerical con- 
trol at very low temperature, compared to the critical temperature of the phase transition, 
we can not plot all points in superconducting phase (the lowest temperature in Figure(6) 
and Figure(8) is about 0.64T C and 0.68T C separately). However, it is not obtrusive to con- 
clude that the entanglement entropy does behave as thermal entropy at least for not very 
low temperature. Furthermore, as the thermal entropy is equivalent to the black hole en- 
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Figure 8: The entanglement entropy in a = 0.447 function of temperature at 

fixed belt width (left plot for £/i = 4 and right plot for £/i = 6). The dashed blue curves are 
from the RN-AdS solutions, while the solid curves are from the superconductor solutions. 
The physical curve is determined by choosing the dashed blue curve above T c , indicated 
by the vertical dotted line, and the curve which has the lowest entropy below T c . 



tropy in the holographic setup, our calculation seems to support the viewpoint that black 
hole entropy is due to the entanglement entropy [40, 41, 42, 43]. 

5 Conclusion and discussions 

In a recent paper [33], we reported the behavior of entanglement entropy in the holo- 
graphical insulator /superconductor phase transition, where a non-monotonic behavior of 
the entanglement entropy was found as the change of chemical potential. Because of the 
absence of horizon in the soliton background, both the temperature and thermal entropy 
in dual boundary system do vanish. It is unable to extract the relationship between en- 
tanglement entropy and thermal entropy or black hole entropy. In this paper we overcome 
the shortcoming by studying the holographic p-wave superconductor at finite temperature. 
The other motivation is to see how the behavior of the entanglement entropy will change as 
we tune the parameter a, the ratio of the gravitational constant to the Yang-Mills coupling 
constant, which can change the order of the phase transition. 

In the fully back reacted case, we found ^ scales as T — T c near critical point, which is 
the same as the result in the probe limit [14]. As noted in Ref. [14] this scaling behavior 
is different from the transition between superconductivity and pseudogap state of high T c 
materials, where the fraction is finite and nonvanishing. The behavior is reminiscent of 
superfluid properties of He II. However, the critical exponent in He II is about 0.67, while 
here it is one. 

We found that the qualitative behavior of entanglement entropy is dramatically different 
for sufficiently small and large a cases as we lower the temperature. For the case a < a c , 
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Figure 9: The thermal entropy St as a function of temperature for a = 0.316 (left plot) 
and a = 0.447 (right plot). The dashed blue curves are from the RN-AdS solutions, while 
the solid curves are from the superconductor solutions. A discontinuous slope at T c in the 
left plot obviously indicates a second-order transition, while the jump in the right plot 
indicates a first-order transition. 

the entanglement entropy is continuous at the critical temperature T c , while it has a jump 
for a > a c case. That is, the behavior of the entanglement entropy shows that the phase 
transition is second order when a < a c , while it is first order as a > a c . 

When the belt width is very large, i.e., 7.4 probes deeply into the bulk, the entanglement 
entropy is extensive as the thermal entropy of the bulk black holes. In the opposite limit, the 
behavior perfectly fits the general form obtained from four- dimensional conformal theories. 
Motivated by Ref. [39], there may exists a crossover function connecting the two limits. We 
compared the behaviors of the entanglement entropy and thermal entropy during the whole 
process of the superconductor phase transition. It shows that they behave qualitatively 
the same. This result is nontrivial since the belt width here is neither too large nor too 
small (note that l\i = 4 and £fi = 6 in Figure. (6) and Figure. (8)). From the minimal 
surface picture, the minimal surface 7^4 can not be very close to "hugging" the horizon 
at these intermediate scales. Due to the lack of numerical control at low temperature, 
in our numerical calculation we are not able to display the behavior of the entanglement 
entropy at sufficiently low temperature, compared to the critical temperature of the phase 
transition. As a result, we can conclude that the entanglement entropy is a good probe to 
the holographic pase transition, and that its behavior can indicate not only the appearance, 
but also the order of the phase transition. 

Here it is quite interesting to compare our results with those from the case with holo- 
graphic s-wave superconductors in Ref. [32]. The model studied there is a 5*0(3) x 
5*0(3) invariant truncation of four- dimensional M = 8 gauged supergravity [44]. Dif- 
ferent from other "top-down" models of holographic superconductors (see, for example, 
Refs. [45, 46, 47]), where the superconductor phase transition is a second order one, the 
phase transition studied in Ref. [32] will be second order or first order, depending on the 
boundary condition of a scalar field in the model. In the second order phase transition 
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case, the entanglement entropy shows the feature of the second order phase transition: the 
entanglement entropy is continuous and its slop has a jump at the critical temperature. On 
the order hand, in the case of first order phase transition, there exists not only a jump for 
the entanglement entropy at the critical temperature, but also a kink shape (swallowtail 
shape) for the entanglement entropy for a given range of strip widths in the supercon- 
ducting phase. The kink persists even for zero temperature solutions. It was argued there 
that the appearance of the kink can be attributed to the existence of a new scale in the 
theory and the entanglement entropy is a good probe of the new scale. The new scale can 
be viewed as a finite correlation length. Quite interestingly, in our p-wave superconductor 
model, the kind does not appear even in the first order phase transition case (see Figure. 7). 
This suggests that the existence of this new scale is not universal in holographic supercon- 
ductor models. Note that in this p-wave model, not only the U(l) symmetry but also the 
rotational symmetry are broken, hence, the superconducting phase is anisotropic. On the 
other hand, only the U(l) symmetry is broken in the s-wave model studied in Ref. [32] and 
the superconducting phase is isotropic. We suspect that the anisotropy might be respon- 
sible for the absence of the kink in the superconducting phase in the first order transition 
case. Clearly it would be quite interesting to investigate whether such kink is common in 
other s-wave superconductor models. 

Note that the zero temperature solution studied in Ref. [32] is an RG flow between 
two AdS spaces. Indeed, the existence of the kink for the entanglement entropy in such 
backgrounds has also been observed in Ref. [48] (see also Ref. [49]) where null energy 
condition is used to constrain the monotonic behavior of c-function along RG flows. It 
was shown there that the kink would emerge for particular choice of the asymptotically 
AdS geometry, which is obviously shown in Figure(6) in Ref. [48]. More specifically, the 
geometry is controlled by a parameter R. A kink will emerge for R < R c , while it will 
disappear for the case R > R c . In Ref. [32], it has been also analyzed in some details 
when the kink will appear. It shows that the presence of the kink is due to the particular 
potential in the model which in turn further determines the geometry structure of the 
background solutions. Let us further stress here that the kink of entanglement entropy 
also appears in the AdS soliton backgrounds [38, 50, 33], and the signal of the kink acts as 
the "confinement /deconfiment" phase transition. 

Finally we would like to mention that the holographic p-wave superconductor model 
has been extended to include the Gauss-Bonnet term in Refs. [51, 52] where the condensate 
is found to become harder as the Gauss-Bonnet coefficient grows up. On the other hand, 
the original entanglement entropy proposal has been generalized to include some higher 
derivative corrections, such as the Gauss-Bonnet term, in the bulk [53, 54]. The entangle- 
ment entropy in the Gauss-Bonnet gravity has been calculated for many configurations and 
discovered to give additional contributions [48, 55, 56]. There might be some new features 
for the entanglement entropy in the p-wave superconductor models with the Gauss-Bonnet 
term. We wish to report on the related work in future. 
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